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PROGRAM TTHEi FAST SINE-COSINE SUBROUTINE 
(FLOATING POINT) 



1. Purpose 

To provide a subroutine to compute both the sine and cosine of a 
radian argument with a small relative error. Speed and precision, rather 
than memory requirements, are the primary objectives of this subroutine. 

2. Restrictions 

2.1 To obtain the sine and cosine to at least three decimal places 

8 
of accuracy, the absolute size of the argument must be less than 10 

radians. See 3.2 below. Can ! t be relocated into XX60 (XX any channel). 

3. Method 

3.1 The original angle Is multiplied by a suitable multiplier 
( k/JX for arguments in radians) to transform the angle to units of half 
right-angles, or eighths of a circle. The integral part (except for the 
low order bit) is then removed, yielding an intermediate argument 0* F <■ 2. 
The final argument is then obtained by using either + |f| or + |2-F| , 
whichever is smaller, with the sign chosen in accordance with the sign of 
the original angle. This final argument is then used in a continued 
fraction expansion for the sine; namely: 

r 



sine x * x <d Q + x 



di + e l 
l g^- 

i d 2 + x + Z 2 

v ' d 3 + x d 



1. 



Calling the value of the fraction above P, and Yl - P 2 « R, the 
value of the sine and cosine are determined from the last 3 bits of the 
integer part of the argument expressed in half right-angle units, as follows: 
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Sin Cos 



R -P 
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3.2 Accuracy: The maximum relative error of the subroutine for either 

sine or cosine in the basic range - -1L- -£ x ^ .ii. radians is less than 

h k 

2 x 10 . For arguments outside of these limits, the absolute accuracy of 

either function is limited by the number of bits representing the fractional 

portion of the radian argument, m particular, the maximum absolute error 

of either sin or cos for radian arguments outside the basic range appears to 

E-38 
be less than 3x2 , where E is the binary exponent of the argument. 

k. Usage 

U.l The argument must be a floating point number in the A and R 

registers. 

U.2 Calling sequence Note 

4 • 

TRA 33-512 SR L + 0030.0 

o 

Normal return 
The transfer may be in either half word. 

2. 



k.3 Upon normal return sine X is in the A and R registers, and also 
in 777U and 7775; cosine X is in 7776 and 7777, all numbers being in 
floating point form. 

h»h Error re turn : None 

h.$ Hie routine destroys the L and V loops 

k*6 The routine is relocatable. 

k.7 Options available: For arguments in degrees, change location 

L + 0124 to +55 it055l-02 66030 ^ 

> 1A5 

L Q - + 0122 to -00 00000 -00 00021} 

This avoids the necessity for the programmer to convert an argument 

from degrees to radians before transferring to the sin/cos routine. 

5. Coding Information 

5.1 The routine occupies words 0000 to OI373 or 96 words. 

5*2 Unused locations. A total of 26 locations are unused, as follows: 
Location (octal) Number 



L Q + 002U - 0025 


2 


00I42 - oohh 


3 


OOI46 - 00h7 


2 


0062 -0061* 


3 


0067 - 0070 


2 


0072 - 0073 


2 


0102 - 0105 


h 


0107 - 0113 


5 


0123 - 0125 


3 




26 



3. 



5.3 Krasable locations. Location L + 0076 is used for temporary 
storage. 

5.4 Constants (all floating point) 

Value Location (octal) 

d o 

*1 

d 2 
d 3 
e l 
e 2 



h/rr 



0.78539 81633 97 


L Q + 0022 


-0.15789 09220 72 


+ 0020 


184. 56259 0042 


+ 0060 


22.19509 18999 


+ 0065 


56.26206 97951 


+ 0126 


12090.47453 99 


4- 0026 


Y ■ 1.27323 95447 4 


+ 0121 


1. 


+ oo4o 


2. 


+ 0100 



^.$ Execution time: 284 ms maximum 

6. Checkout 

The routine has been checked over a wide range of values, using 
both degrees and radians as arguments. Comparisons with standard tables 
revealed no errors in excess of those indicated in 3.2. 
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0000.0 
0001 .0 
0002.0 

0003.0 

0004.0 

0005.0 
0006.0 
0007.0 



+j0.777)4.0 FCA 
+35.T7.76.0 FST 



+30.7772. 
+57.7770. 



FCA 
TRA 



+34.7772.O FCS 
+35.7776.0 FST 



+30.7774. 
+57.7770. 

+34.777^ 
+35.7776. 

+34.7772. 
+57.7770. 



FCA 

TRA 

FCS 

FST 

FCS 

TRA 



+30.7772.0 FCA 

+35.7776.0 FST 

+3^.777^.0 FCS 

+57.7770.0 TRA 



0020.0 
0021 .0 
0022.0 
0023.0 
0024.0 
0025.0 
0026.0 
0027.0 



-5O.327O.I 
-41 .5242.0 

-00.0000.0 
-00.0001 .0 

+62.2077.0 
+52.4207.0 

+00.0000.0 
-00.0000.0 

-00.0000.0 
-00.0000.0 

-00.0000.0 
-00.0000.0 

+57.1647.1 
-27.5561.1 

+00.0000.0 

-00.0007.0 



0010.0 +30.0126.0 FCA 

+O5.7776.O FDV 

001 1.0 +04.0020.0 FAD 

+O7.7772.O FMP 

0012.0 +04.0022.0 FAD 

+07.7774.0 FMP 

0013.0 +35.7772.O FST 

+O7.7772.O FMP 

001 4.0 +33.7767.O EXT 

+04.0040.0 FAD 

0015.0 +35.7763.0 FST 

+44.7763.0 FSQ 

OOI0.O +35.7774.0 FST 

+57.7771.0 TRA 

0017.0 -77.7777.' 

+77.7777.1 — 



OO3O.O 



0031.0 
0032.0 
OO33.O 
0034.0 
OO35.O 
OO36.O 
OO37.O 



+I5.OO76.O SAX 

+57.0045.0 TRA 

+34.7774.O FCS 

+06.7772.1 FSB 

+04.7774.0 FAD 

+04.7772.1 FAD 

+5O.7766.I TZE 

+35.777^.0 FST 

+43.OOOO.O XAR 

+5O.7767.O TZE 

+5I.7767.O TMI 

+53.OOOO.O XAR 

+06.0100.0 FSB 

+35«777+-0 FST 

+64.0130.0 CTL 

+57. 7760. TRA 



0040.0 
OOUl .0 

oo42.o 
0043.0 
oo44.o 
0045.0 
0046.0 
0047.0 



+ to. 0000.0 

-00.0000.0 

+00.0000.0 
-00.0000.1 

-00.0000.0 
-00.0000.0 

-00.0000.0 
-00.0000.0 

-00.0000.0 
-00.0000.0 

+64. 0050.0 
+57.77fo.o 

-00.0000.0 
-00.0000.0 

-00.0000.0 
-00.0000.0 



CTL 

TRA 



0060.0 
006). 
0062.0 
0063.0 
0064.0 
0065.0 
0066.0 
0067.0 



+56.1100.0 
-27.1*651 ..1 

+00.0000.0 
-00.0004.0 

-00.0000.0 
-00.0000.0 

-00.0000.0 
-00.0000.0 

-00.0000.0 
-00.0000.0 

+54. 3076.0 
-61.2744. 

+00.0000.0 
-00,0002.1 

-00.0000.0 
-00.0000.0 



0050.0 
0051.0 
0052.0 
0053.0 

0054.0 
0055.0 
0056.0 
0057.0 



+66.OO7O.O CTV 

+01.7771.0 ADD 

+6O.777O.O STO 

+OO.7776.O CLA 

+07.0121.0 FMP 

+35.7772.0 FST 

+51.7764.1 TMI 

+04.0115.0 FAD 

+57.7765.O TRA 

+04.0117-0 FAD 

+33-7771.0 EXT 

+01.0106.0 ADD 

+33.01 iH.O EXT 

+60. 777 1.0 STO 

+64. OO3O.O CTL 

+57. 7761.0 TRA 



0070.0 -00.0000.0 — 

-00.0000.0 — 

0071 .0 +35.7774.0 FST 

+57.0000.1 TRA 

0072.0 -00.0000.0 — 

-00.0000.0 — 

0073.0 -00.0000.0 — 

-00.0000.0 — 

0074.0 +4o.oooo.o — 

-00.0000.0 — 

0075.0 +00.0000.0 — 

-00.0023.1 — 

0076.0 -00.0000.0 — 

-00.0000.0 — 

0077.0 +00.0000.0 

-00.0007.0 — 



0100.0 +4o. 0000.0 

-00.0000.0 — 

0101 .0 +00.0000.0 
-00.0001.0 

0102.0 -00.0000.0 — 

-00.0000.0 — 

0103.0 -00.0000.0 — 

-00.0000.0 -— 

0104.0 -00.0000.0 — 

-00.0000.0 

0105.0 -00.0000.0 — 

-00.0000.0 — 

0106,0+64.0000.0 CTL 

+57.7761 .0 TRA 

0107..0 -00.0000.0 — 

-00.0000.0 — 



0120.0 
0121.0 
0122.0 

0123.0 

0124.0 

0125.0 
0126.0 
0127.0 



+00.0000.0 

-00.0023.0 

+50. 5746.0 
-15.5625.0 

+00.0000.0 
-00.0000.1 

-00.0000.0 
-00.0000.0 

-00.0000.0 
- 00 o 0000.0 

-00.0000.0 
-00.0000.0 

+70.2061.0 
+60.0616.0 

+00.0000.0 

-00.0003.0 



0.1 10.0 -00.0000.0 — 

-00.0000.0 — 

0111.0 -00.0000.0 — 

-00.0000.0 — 

0112.0 -00.0000.0 ,— 

-00.0000.0 — 

0113.0 -00.0000.0 — 

" -00.0000.0 — 

0114.0 +77.77TZ.1 — 

+77.7766.0 — 

0115.0 +40o 0000.0 — 

-00.0000.0 — 

0116.0 +00.0000.0 — 

-00.0023.0 — 

0117.0 +60.0000.0 

-00.0007.0 — 



0130.0 
0131.0 
0132.0 

0133.0 
0134.0 

0135.0 

0136.0 

0137.0 



+30.7772. 

+52.7762. 

+34.7774. 

+60. 7774. 

+30.7774. 
+07.7774. 

+35.7772.. 

+04.0065. 

+35.776: 

+30.002* 

+05.7763. 
+04.7772. 

+04.0060. 

+35.7776. 

+64.0010, 

+57.77600 



FCA 

TPL 

FCS 

STO 

FCA 

FMP 

FST 

FAD 

,0 FST 

FCA 

FDV 

0* FAD 

FAD 

0' FST 

CTL 

0" TRA 



